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The merging procedure of tree-level matrix elements and the subsequent parton shower as imple¬ 
mented in the new event generator SHERPA will be validated for the example of single gauge boson 
production at the LHC. The validation includes consistency checks and comparisons to results ob¬ 
tained from other event generators. In particular, comparisons with full next-to-leading order QCD 
calculations prove SHERPA’s ability to correctly account for additional hard QCD radiation present 
in these processes. 

PACS numbers: 13.85.-t, 13.85.Qk, 13.87.-a 


I. INTRODUCTION 

The production of electroweak gauge bosons, which decay leptonically, is one of the most prominent examples for 
hard processes at hadron colliders and one of the first applications of perturbative (^D in such reactions. In fact, 
the next-to-leading order (NLO) corrections to this process in QCD, calculated by Ulliill provided the first 
calculation of such corrections for hadron collisions. Later, their production cross section has been calculated at 
NNLO by 1^0. Recently, the first distribution determined at NNLO related to these processes, namely the boson 
rapidity, has been calculated by Q. In addition, there is a large number of computer programs dealing with single 
gauge boson production. They range from RESBOS |^, which resums soft gluon effects in these processes, to codes, 
that evaluate cross sections at the LO level for t he p roduction of gauge bosons accompanied by jets. Examples for the 
latter include specialised ones, such as VECBOS [^, and general ones, usually called parton level generators, such as 
COMPHEP GRACE/GROPPA [il[l3|, MADGRAPH/MADEVENT [ll[l3, ALPGEN [Tg, and AMEGIC++ [13. Furthermore, the 
first package called MCFM has been made available that calculates total and differential cross sections at NLO precision 
for the production of gauge bosons with up to two jets 0,^|. This reflects the importance of this particular process. 
At the LHC becoming operational in the near future, the gauge bosons will be produced with unprecedented rates. 
For instance, at luminosities of £ « lO^^cm^/s the production and leptonic decay of a single W boson will occur with 
a frequency of around 20 Hz, rendering this process a prime candidate for luminosity monitoring miilllllil. Of 
course, these large rates will allow to measure the gauge bosons parameters, such as their masses and widths, with a 
precision [23,123 beyond what could be reached at previous collider experiments [t^ 1^ 1^ 1^ 1^ 1^ 1^ 1^ 1^ . 
At the CERN LHC, the production of W and Z bosons together with jets also constitute an important background 
to all kinds of searches for new physics, such as supersymmetry. An example for this is the production and decay of 
gluinos, where the production of jets plus a Z boson decaying into neutrinos forms a major background. 

In a previous analysis [13 it has been shown that some results for this type of process, i.e. the production of single 
gauge b osons plus extra jets, as obtained by other multi-purpose event generators such as PYTHIA [33[13i HERWIG 
[SI l4l| or even MCONLO differ significantly from the results obtained by SHERPA ^3- In particular, it has 

been shown that already at the Fermilab Tevatron, operating at roughly 2 TeV centre-of-mass energy, the additional 
jets are produced at significantly larger transverse momenta. The reason for this difference is the way the different 
codes implement the knowledge of exact matrix elements for the production of multi-particle final states. In both, 
PYTHIA and HERWIG, they are included at first order in as through a correction of the first hard emission on the 
corresponding qq Vg oi qg ^ Vq matrix element, where V stands for the vector boson [dR ItR Id^ . In MCSNLD 
the full first order correction, including both, the virtual and the real parts, are matched with the parton shower. 
This has the additional benefit that MCSNLO reproduces correctly the total production rate of single gauge bosons and 
the spectrum of the first additional jet at first order in as- In contrast, in SHERPA a method has been implemented 
that consistently adds different matrix elements at the tree level for different jet multiplicities and merges them 
with the parton shower. The basic idea in this approach is to internally define a region of jet production (hard 
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parton emissions) and a region of jet evolution (soft parton emissions). The two regimes are divided by a A:_L-type 
of jet measure [dol 1^ . Leading higher order effects are incorporated by reweighting the matrix elements with 
appropriate Sudakov form factors. Formal independence at leading logarithmic order of the overall result on the 
jet measure is achieved by suitable starting conditions and vetoing hard emissions inside the parton shower. This 
approach was presented for the first time for e+e“ collisions; it has been extended to hadronic collisions in [h^. 
A reformulation for dipole cascades has been presented in [^ . The algorithm is implemented in a fully automated 
way and in full generality in SHERPA, some other realisations 1^ proved the flexibility and the validity of the 
approach. 

In this publication, the previous analysis will be extended to the case of the CERN LHC, operating in the pp 
mode at 14 TeV centre-of-mass energy. In Sec. QH a number of consistency checks will focus on the independence of 
the results on variations of the internal jet definition and of the number of matrix elements involved. Also, the effect 
of scale variations in both the matrix elements and the parton shower is investigated there. Then, in Sec. EH results 
obtained with SHERPA will be contrasted to those obtained from fixed order (LO and NLO) calculations provided by 
MCFM. Following this, different multi-purpose event generators, namely PYTHIA, MCQNLO and SHERPA, will be compared 
in Sec. im before the conclusions will summarise the findings. 


II. CONSISTENCY CHECKS 

Before comparing the results of SHERPA with those of other programs, some consistency checks will be performed. To 
do so, the dependence of some observables in reactions of the type pp —*■ e“'"e“ -I- A on internal parameters intrinsic 
for the merging procedure will be investigated. In particular, these parameters are the internal jet resolution cut 
Qcut and the maximal number Umax of final state partons (giving rise to jets) described through matrix elements. 
The former parameter defines the transition of the matrix element domain to the phase space region covered by 
the parton shower during event generation. In principle, the actual value of this parameter can be chosen freely, 
nevertheless their exist criteria that guide such a choice. For very low values of Qcut the evaluation of the matrix 
elements becomes very challenging and potentially inefficient once jet cuts are performed on the analysis level harder 
than the generation cut Qcut- The upper limit is defined by the scale where jets produced by the parton shower 
start to disagree significantly from such produced by equivalent matrix elements. To study especially the effect of 
the upper limit, the values used in this analysis will range over nearly one order of magnitude, from 15 to 100 GeV. 
The choice of the number of matrix element partons taken into account may be steered by two aspects. First of all, 
JT-max has to be sufficiently large to properly account for the phase space region the observable under consideration is 
sensitive to. As an example, consider the transverse momentum of the boson compared to that of the, say, third jet. 
It is obvious that a rather inclusive quantity such as the former may be appropriately described with lower values of 
?^max than the latter observable. On the other hand, the upper limit on Timax is given by the availability of the matrix 
elements at all and by the potentially large amount of CPU time the evaluation of multi-leg matrix elements requires. 
Within SHERPA matrix elements with up to four extra partons can be delivered for the processes under consideration 
in this publication. After evaluating the sensitivity of the results on the principal parameters defining the merging 
procedure, Qcut and Umax, the effect of scale variations will be investigated. This, together with the dependence on 
Qcut and rimax yields an estimate for the uncertainty related to predictions of SHERPA. 

The results presented in this section were generated with the following setups: when varying the jet resolution 
parameter Qcut, the maximal number of final state partons Umax has been set to rtmax = 3. When studying the impact 
of different matrix element multiplicities, the scale Qcut has been fixed to Qcut = 15 GeV; this clearly maximises the 
impact of the higher order matrix elements. When scale variations are under consideration, the choice Qcut = 20 GeV 
and n„iax = 2 has been made. 

In the following, .Z-boson production will be investigated in more detail. Nevertheless, the process under consideration 
is pp —!■ Zj^* e~ A, where the full ^-Z interference is taken into account and spin correlations are fully 

respected. Further input parameters used and the phase space cuts applied are summarised in the appendix^! Note 
that the cut on the invariant mass of the lepton pair is just ruee > 15 GeV which is rather small. The description 
of such low mass lepton pairs constitutes a real challenge for the description through the merging prescription. The 
reason is that at large Qcut = 0(100 GeV), lepton pairs with such low invariant mass clearly are softer than any jet 
produced through the matrix element, rendering a consistent merging a complicated task. 


A. Observables related to the leptons 

Starting from more inclusive observables, first of all the effect of parameter variation on lepton observables will be 
considered. In Fig.E the Pa. spectra of both the lepton pair (upper row) and of the electron alone (lower row) are 
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FIG. 1: p_l{Z) (upper row) and p±(e ) (lower row) for Qcut = 15 GeV, 50 GeV and 100 GeV (from left to right). The dashed 
reference spectrum has been obtained after averaging the results for Qcut = 15, 20, 30, 50, 100 GeV. 


shown for three different values of Qcut: from left to right, in the columns Qcut = 15, 50, 100 GeV, as indicated 
by the thin vertical lines. In each plot, the resulting spectrum is compared to a reference obtained from averaging 
the results for Qcut = 15, 20, 30, 50, 100 GeV. In this and all other plots, contributions stemming from the different 
matrix element multiplicities are indicated through coloured lines. 

Using Qcut = 15 GeV obviously produces the hardest boson/lepton spectrum. It is the smallest cut considered here 
and therefore the distributions are dominated by matrix elements that in contrast to the parton shower favour rather 
hard parton kinematics. For very high pj_ the distributions are almost completely covered by the matrix element with 
the highest multiplicity (n-max = 3)- This shows that the LHC provides enough phase space to produce a sufficient 
amount of events with three and more jets of p_L > Qcut- For the case of Qcut = 50 GeV the situation is slightly 
different. The high-pj_ tail is filled to an equal amount by the different multiplicities, the total sum being slightly 
below the reference curve. This reference curve contains three results with jet resolutions smaller than 50 GeV that 
somehow dominate the averaged result. The spectrum for Qcut = 100 GeV starts to underestimate w.r.t. the reference 
the boson transverse momentum at pj_ « 35 GeV and the lepton pj_ for values larger than 60 GeV. To understand this, 
one has to remember that the boson transverse momentum for values below the resolution cut is almost completely 
covered by the parton shower. The shower description, however, is known to suffer from a lack of hard QCD radiation. 
This leaves not enough hard partons, the boson can recoil against. Beyond this influence of the Qcut variation on 
the intermediate and high boson transverse momenta, it has to be noted that all curves are very smooth around the 
jet resolution cut. Although the cut defines a rather sharp transition from the parton shower to the matrix element 
domain no signihcant holes in the boson and lepton p±^ spectra can be observed. 

In Fig. [21 the pseudo-rapidity spectra of the lepton pair and the single electron are displayed; again for Qcut = 
15, 50, 100 GeV with the same way of generating the reference. While the electron observable is nearly unaltered, the 
differences in the ry distribution of the lepton pair can be understood easily: the smaller the chosen cut, the larger the 
influence of the matrix elements with extra external legs. These matrix elements however prefer to produce the boson 
much more central than the parton shower does. This effect yields slightly tighter spectra with the central rapidities 
being pronounced for smaller resolution cuts. 

The effect of varying Umax on the transverse momentum and pseudo-rapidity spectrum of the lepton pair is exhibited 
in Fig. 0 In this figure, results are compared for rimax = 2, 3, 4. In each plot, a reference result is given with the 
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FIG. 2: T]{Z) (upper row) and T]{e ) (lower row) for Qcut = 15 GeV, 50 GeV and 100 GeV (from left to right). The dashed 
reference spectrum has been obtained after averaging the results for Qcut = 15, 20, 30, 50, 100 GeV. 


corresponding = n„iax — 1. For the case of the p± distribution it has already been observed that the high-p^ 
region is described through higher multiplicity matrix elements. As a consequence, it is the high-p_L region that is 
affected by the variation of Umax- However, while the effect is clearly noticeable when going from one to two extra 
partons the change becomes smaller the more matrix elements are included. From the very right plot one can conclude 
that considering Z + S extra parton matrix elements is a reasonable choice to simulate inclusive Z production. The 
change in the ry distribution for different rimax is as expected, considering what has already been seen for varying the 
jet resolution. The higher multiplicity matrix elements favour the region of small |ry| yielding slightly tighter pseudo¬ 
rapidity distributions. Again, the more matrix elements have been taken into account the smaller the influence when 
adding an even higher multiplicity. 

In comparison to what has been observed when studying gauge boson production at the Fermilab Tevatron , the 
LHC provides much more phase space for additional hard QCD radiation, enhancing the influence of higher order 
matrix elements. Therefore a modest value of the jet resolution parameter and the inclusion of a sufficient large 
number of matrix element legs is advisable for LHC analyses. 


B. Jet observables 

As has already been seen in the previous publication [^, a very sensitive test of the merging procedure is provided 
by observables based on jets. In particular, differential jet rates have turned out to be very useful, since they clearly 
show how the matrix elements and the parton showers interact in filling the phase space below and above the jet 
resolution cut. In Fig. ^ differential jet rates using the Run H k± clustering algorithm with i? = 1 are depicted. 
They signal the relevant Q value of the fcj_-algorithm, where an (n -|- l)-jet event turns into an n-jet event. Again, 
the results for three different values of Qcut are depicted: from left to right, in the columns Qcut = 15, 30, 100 GeV, 
as indicated by the thin vertical lines. In each plot, the resulting spectrum is compared to the average of the results 
for Qcut = 15, 20, 30, 50, 100 GeV. In the three rows, the differential jet rates for the 1^0, the 2^1, and for the 
3 —t 2 transition (from top to bottom) are shown. Starting the discussion with the results for Qcut = 30 GeV, very 
good agreement with the reference curves can be observed. While the 3 —> 2 transition is very smooth around the cut 
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FIG. 3: p±{Z) (upper row) and ^{Z) (lower row) for Qcut = 15 GeV and different maximal numbers (2-4, from left to right) 
of ME jets included. The dashed line corresponds to the maximal number of ME jets reduced by one. 


the results for 1 ^ 0 and 2^1 exhibit small dips at the cut scale. Since the kinematics of the matrix elements is 
altered when the parton shower is attached, mismatches of the parton configurations close to the cut occur, leading 
to the dips. Similar structures can be observed for the case of Qcut = 100 GeV. However, more obvious here is that 
the parton shower fails to fill the phase space for hard emissions up to this very large cut. For Qcut = 15 GeV no 
visible dips at the cut scale are observed. Instead, this sample seems to slightly overestimate the contributions from 
higher order matrix elements w.r.t the reference. A small kink at Qcut can be observed for the 1^0 and 2^1 
transition. These residual dependences of the results on Qcut may be used to tune the perturbative part of the Monte 
Carlo event generator. 

In Fig. 0 the p± spectrum of the jet in exclusive Z -b Ijet production is shown for three choices of the jet resolution 
scale, Qcut = 15, 50, 100 GeV, indicated by the thin vertical line. The results are contrasted with a reference curve, 
again the average of results for Qcut = 15, 20, 30, 50, 100 GeV. The jet has been defined using the Run II fc_L-algorithm 
with a minimal jet-p_L of 20 GeV and R = 0.4. The smallest value of Qcut presented here, namely 15 GeV, is smaller 
than the actual jet cut used in the analysis. Accordingly, matrix elements with more than one extra leg have a 
non-vanishing influence on the jet-pj_ distribution. This changes as soon as Qcut becomes larger than 20 GeV. For 
Qcut = 50 GeV and even more for Qcut = 100 GeV the contributions from matrix elements with rimax > 1 are almost 
negligible. There, only a small dip in the pj_ distribution around the resolution scale can be observed. As has been 
seen in the transverse momentum distribution of the lepton pair, cf. Fig.^ for Qcut = 100 GeV, the shower is not able 
to fill the full phase space below the cut properly. However, the overall agreement of the three results is satisfactory, 
keeping in mind the large parameter range used for Qcut- 

To highlight the effect of taking into account different maximal numbers of final state partons through matrix elements, 
a two-jet correlation is exhibited in Fig. El There, the relative transverse angle A(/) between the two hardest jets in 
inclusive Z + 2jet production is displayed; from left to right, rimax has set to Umax = 2, 3, 4. Each result is contrasted 
with a reference that has been obtained with = rimax ~ 1- From the very left plot it is clear, that the one- 

jet matrix element is incapable of correctly describing the A;/) distribution since the parton shower does not treat 
interferences properly. On the other hand, as soon as Umax > 2, the two-jet correlations are consistently described 
and changes due to the inclusion of higher order matrix elements are rather modest. 
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FIG. 4: Differential jet rates for the 1 ^ 0, 2 —> 1 and 3^2 transition (top to bottom), for Qcut = 15 GeV, 30 GeV, and 
100 GeV (from left to right). The dashed reference curve in each plot is obtained after averaging the corresponding results for 
= 15, 20, 30, 50, 100 GeV. 


C. Variation of renormalisation and factorisation scale 

The algorithm as implemented in SHERPA determines the renormalisation and factorisation scales used in a specific 
calculation. Of course there is some intrinsic freedom in defining the scales used for the evaluation of the PDF or 
the strong coupling constant. In particular, the scales used can be multiplied by constant factors, as long as this 
alteration is applied both in the matrix element evaluation and reweighting and in the parton shower. This restriction 
is due to the construction of how the leading logarithmic dependence on Qcut is eliminated. The dependence of the 
SHERPA results with respect to such scale variations is studied in Fig. 0 and Fig.|Sl Results obtained with the default 
scale choices are confronted with results obtained when all scales appearing in the coupling constants and PDFs are 
multiplied by common factors of 0.5 and 2. 

In Fig.[7|the transverse momentum and pseudo-rapidity distribution of the Zj^* boson are depicted. For the case of 
the p_L spectrum, except for the very first bins the spectrum obtained with a factor of 0.5 (2) is always above (below) 
the default result. The differences are rather constant and of the order of 10— 15%. As has been seen before, cf. Fig.^ 
and 13 for transverse momenta above the cut scale the distribution is predominantly described by hi gher order matrix 
elements, whose scale dependence is known to be tamed with respect to the lowest order process Il0l| . This lowest 
order process, however, dominates the region of very low boson momenta. There the 2^2 cross section exhibits a 
strong decline when the scales become smaller. This effect potentially leads to the reversal of the discrepancies in 
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FIG. 5: p± of the jet in exclusive Z + Ijet production. For the jet definition, the Run II fcx-algorithm with R = 0.4 and 
> 20 GeV is used. From left to right, results for Qcut = 15, 50, 100 GeV are contrasted with a reference: the average of 
the results for Qcut = 15, 20, 30, 50, 100 GeV. The thin horizontal line indicates the jet resolution scale used. 



FIG. 6: A(^i 2 for Qcut = 15 GeV and different maximal numbers of ME jets included. The dashed line corresponds to the 
reference result obtained with njjlx = fimax — 1. 

the soft region. The situation in the case of the pseudo-rapidity distribution is very similar. From Fig. Inland Fig. O 
one can read off that the region of large values of |? 7 | is described by the parton shower attached to the 2^2 matrix 
element. For |? 7 | > 5 the spectrum, where all scales have been multiplied by a factor of two, is enhanced up to 20%. 
A factor of 0.5 on the other hand depopulates this phase space region by up to 20%. In the intermediate range of 
pseudo-rapidity the deviations of the two spectra from the default scale choice are well below 10%. 

In Fig. IHl the transverse momentum distribution of the hardest jet in inclusive Z production is depicted. In 
contrast to the two distributions above, this result has no significant contribution from the leading order 2 —> 2 
process. Therefore, the two results obtained after scale manipulation do not cross each other. Over the whole range 
of jet transverse momentum the deviations of the two curves from the default result are very moderate. 

It can be concluded that the predictions of SHERPA show rather mild variations over a wide range of the phase space 
when multiplying all scales appearing in the coupling constants and PDFs by common factors of 0.5 and 2. The 
largest deviations from the default choice of scales are observed in those phase space regions that are predominantly 
covered by the 2—^2 matrix element with the parton shower attached. 


III. SHERPA VS. NLO RESULTS 

Having investigated the self-consistency of the merging procedure as implemented in SHERPA its parton level results 
are compared with those from MCFM, v. 4.0, For the class of processes studied here, MCFM is capable 

to calculate total and fully differential cross sections at next-to-leading order in the strong coupling constant for 
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FIG. 7: The transverse momentum (left) and pseudo-rapidity (right) distribution of the boson in inclusive production at 

the LHC and their dependence on different choices for the factorisation and renormalisation scale. Results have been obtained 
using Qcut = 20 GeV and rimax = 2. The black solid lines indicate the default hadron level result of SHERPA. To obtain the 
green dashed (red dotted) line, all scales appearing in the coupling constants and the PDFs in both the matrix elements and 
the parton showers are multiplied by a factor of 0.5 (2). In the lower parts of the plots the variations of the results with respect 
to the default scale choice are presented. 



p.p (first jet) [GeV] 


FIG. 8: The px spectrum of the hardest jet in inclusive Zj^* production at the LHC and its dependence on different choices for 
the factorisation and renormalisation scale. Results have been obtained using Qcut = 20 GeV and Umax = 2. Jets are defined 
through the Run H fex-algorithm with R — 0.4 and > 20 GeV. The black solid line corresponds to the default hadron level 
result of SHERPA. To obtain the green dashed (red dotted) line, all scales appearing in the coupling constants and the PDFs 
in both the matrix elements and the parton showers are multiplied by a factor of 0.5 (2). In the lower part of the plot the 
variations of the results with respect to the default scale choice are presented. 


(Z/ 7 * ^ l^l~) -b 0,1, 2 and (W^ lui) -t- 0,1, 2 partons. For all calculations with MCFM the cteq61 PDF has been 
used, and as(jnz) = 0.118 in accordance with the value of the PDF evolution. The renormalisation and factorisation 
scales have been chosen to be identical with the bosons mass, i.e. = rnz or ub = fip = rnw, respectively. 

Phase space cuts are listed in the appendix IXI In contrast to a previous publication l37l . this time only “inclusive” 
quantities are compared. For the next-to-leading order calculation this translates into an unconstrained phase space 
for the real higher order correction. Thus, the higher order corrections may give rise to an additional jet. The SHERPA 
results were obtained after the parton shower evolution. For the sake of a better comparison, all curves have been 
normalised to one, eliminating the enhancement of the cross section due to the NLO corrections. 

First of all, in Figs. M and cni the p± spectra of the hardest jet in inclusive FF+ + Ijet, W + Ijet, and Zj^* Ijet 
production are exhibited. For all cases, results at leading and at next-to-leading order were contrasted with results 
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FIG. 9: The pr distribution of the hardest jet for inclusive W (left panel) and (right panel) plus one jet events at the 
LHC. 


Zj + X @ LHC 



FIG. 10: The pr distribution of the hardest jet for inclusive Zj'^* plus one jet events at the LHC. 


from SHERPA that were obtained with rimax = 1 and with rtmax = 2, respectively. In all plots the high-p_L tail is 
significantly enhanced when going from LO to NLO. The SHERPA samples with rimax = 2 show the same behavior 
but tend to pronounce the high-pj_ region even more. This is in striking contrast to the rimax = 1 samples. They 
are incapable of recovering the shape of the distribution at NLO, and tend to look like the LO result. This is not 
surprising. The NLO calculation takes into account tree-level matrix elements with two final state partons as the real 
contribution to the NLO result. Due to the large phase space available at the LHC this real contribution tends to 
produce an extra jet that alters the kinematics of the first jet. Obviously this significant change in the kinematics can 
not be appropriately recovered by the parton shower. The Umax = 2 SHERPA samples also include the parton shower, 
resulting in increased parton emission thus enhancing the high-p^ tail even more. It would for sure be instructive to 
check this behaviour with a resummed NLO computation for these processes. 

In Figs. E] and El the p± spectra of the two hardest jets in inclusive VF+ + 2jet, W~ + 2jet, and Zj^* -P 2jet 
production are displayed. This time next-to-leading order results from MCFM are compared with the corresponding 
SHERPA samples with Umax = 2. It has been shown in M that the shapes of the distributions when going from LO 
to NLO are quite stable. The slopes of the next-to-leading order and the SHERPA result are in good agreement, the 
latter having the tendency to produce the first jet slightly harder. In Fig.Elthe p± spectra of the two hardest jets in 
inclusive W~ + 2jet production are displayed once more. This time, however, the renormalisation and factorisation 
scale in the NLO calculations has been chosen as pa = pp = 2 m-yv. For this choice of the scales the agreement of MCFM 
and SHERPA is even better. This highlights the effect of scale variations, a good way to estimate residual uncertainties 
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FIG. 11: The pr distribution of the two hardest jets in inclusive W (left panel) and (right panel) plus two jet production 
at the LHC. 


Zjj + X @ LHC 



FIG. 12: The pr distribution of the two hardest jets in inclusive ^/7* plus two jet production at the LHC. 


due to higher order corrections and shows that the results of SHERPA are well within theoretical uncertainties 


IV. SHERPA VS. MC@NLO AND PYTHIA 

In this section^ hadron-level results of SHERPA will be compared with those of two other event generators, namely 
MCQNLO and PYTHIA [^. The former program incorporates a consistent matching of a full-fledged 

next-to-leading order calculation with the parton shower provided by HERWIG ^^^ns employs an angular- 

ordered shower, taking full account of coherence effects. In contrast, PYTHIA uses tree-level matrix elements, in this 
case for qq e“'"e“ and it employs a virtuality-ordered parton shower to model further emissions. In this framework, 
coherence effects are approximated through an explicit veto on rising opening angles in the splitting. Hence, the 
parton shower implementations of PYTHIA and SHERPA are quite similar. However, in order to account in PYTHIA for 
jets with a p± larger than the “natural” starting scale of the parton shower equal to the invariant mass of the lepton 


^ It should be noted that the effect of this scale variation on the total cross section merely is of the order of 1%, although the shape of 
the distribution in the high-py tail changes considerably. 
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FIG. 13: The pt distribution of the first and second hardest jet in inclusive W plus two jet production at the LHC. For the 
NLO calculation the renormalisation and factorisation scale have been chosen to fiR = /if ~ 2mw. 


pair, the starting scale has been increased to the centre-of-mass energy of the proton-proton system, i.e. to 14 TeV. 
This choice is supplemented with a matrix element correction procedure implemented through reweighting meant to 
reproduce the exact matrix element for the emission of an additional jet. The precise setups for both codes can be 
found in appendix 0 

First of all, the results of the three programs for some rather inclusive quantities are compared. The transverse 
momentum and pseudo-rapidity distributions of the produced bosons are presented in Fig. 1141 and Fig. 1151 The 
SHERPA predictions depicted have been obtained with rimax = 1 and Qcut = 20 GeV, in order to match the approaches 
of the other codes. In order to compare the different samples, they all have been subject to a cut on the boson 
invariant mass of the form 


my — 30 • Fy < my < my + 30 • Fy , (1) 

where no additional phase space cuts have been applied. All distributions have been normalised to their respective 
cross section. 

The results for both processes look very similar. The boson transverse momentum distributions of MCONLO and SHERPA 
agree fairly well. In the case of Zj^* production they match nearly perfectly for values of p± > 100 GeV. In the 
intermediate range of 10 GeV < p± < 100 GeV SHERPA apparently is below MCSNLO. This discrepancy may have its 
origin in the different shower approaches used within the two programs. This statement is also hinted at by the fact 
that the PYTHIA result follows the SHERPA distribution for p^ < 35 GeV. For larger values of p±, however, the PYTHIA 
distribution is far below MCSNLD and SHERPA predicting much less bosons with large transverse momentum. For the 
case of W~^ production the MCSNLD and SHERPA prediction cross at p± ~ 60 GeV. SHERPA produces slightly less events 
with smaller boson pj_ and tends to pronounce the high p± region a bit. Again PYTHIA produces fewer bosons with 
intermediate and large boson transverse momenta. Looking at the pseudo-rapidity distributions, it can be recognized, 
that MCONLO and SHERPA both tend to produce the bosons much more central than PYTHIA. Especially the region of 
\ri\ < 4 is filled significantly with respect to PYTHIA, which, in contrast, features a much broader shape. This effect 
is of course directly correlated to the larger amount of hard QGD radiation the other two programs produce, since 
this enhanced QGD radiation allows for larger boson recoils. Moreover, from Fig. 0 it can be anticipated how the 
SHERPA results change under the inclusion of matrix elements with extra QGD legs: the boson transverse momentum 
distribution develops a more pronounced large p± tail and the very cental region of ry is filled even more, thus reducing 
the amount of events with large values of [Tyj. So while the p± spectra would be slightly harder than those of MCSNLD 
the 77 distributions would fit even better than they do for the case of including V -I- 0 and V + 1 parton matrix elements 
only. 

For the comparison of jet observables, only the case of Z production is studied. The qualitative statements implied 
by it, however, will hold true as well in the case of W production. To judge the abilities of the three programs to produce 
extra hard QGD radiation associated to the electro-weak gauge bosons, the transverse momentum distributions of 
the three hardest jets are depicted in Fig. Hi and Fig. El For this comparison in addition to the cut on the boson 
transverse mass according to Eq. o the jet criteria and phase space cuts of appendix IXlhave been applied. For SHERPA 
the jet resolution parameter has been set to Qcut = 20 GeV. The standard sample for this comparison uses again only 
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FIG. 14: The p± (left) and rj (right) distribution of the lepton pair in inclusive production of a Z/ 7 * boson decaying into e'’'e 
at the LHC. The results of the generators MCSNLO (red), PYTHIA (blue), and SHERPA (black) are compared. 
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FIG. 15: The W'^ p± (left) and p (right) distribution in inclusive production at the LHC. The results of the generators MCSNLO 
(red), PYTHIA (blue), and SHERPA (black) are compared. 


matrix elements with up to one additional parton. To test the predictions of SHERPA samples with Umax = 2(3) have 
been considered as well, the corresponding results are shown as dashed (dotted) lines in the plots. Since it is actually 
the production rate that is important here, this time the curves have not been normalized. Instead the corresponding 
differential cross sections are presented. 

For the hardest jet the predictions of MCSNLO and SHERPA agree rather well. The total rate of SHERPA is 12% smaller 
than that predicted by MCSNLO, the distribution, however, has a slightly harder tail. This difference in rate can be 
traced back to the different inclusive production cross sections. However, for Umax = 2, the two total cross sections of 
Z + Ijet nearly coincide (cf. the dashed black curve in Fig. lltill . In terms of shape, SHERPA apparently favours jets with 
larger transverse momentum. As has been seen in the closely related boson p± distribution in Fig.^^ PYTHIA predicts 
a much smaller rate (60% w.r.t. the rate predicted by MCQNLO) for the production of extra hard QCD radiation with 
a softer distribution. 

For the second jet the situation changes significantly. Here, even in the case of including only matrix elements with 
up to one extra parton the two-jet rate predicted by SHERPA is 17% larger than that of MCSNLD. Including matrix 
elements with two extra partons the difference becomes nearly 90%. As for the case of the first jet, PYTHIA predicts 
the radiation of a second jet with a much smaller rate. Similar statements hold true when looking at the third jet but 
this time the differences are even larger. Note, that a reliable prediction of the three jet rate requires the inclusion 
of matrix elements with at least three extra partons. While the sample with matrix elements up to one extra parton 
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FIG. 16: The px distribution of the first (left) and second (right) hardest jet in inclusive production at the LHC as 

obtained by MCSNLO (red), PYTHIA (blue) and SHERPA (black). The dashed curves correspond to the predictions of SHERPA when 
matrix elements for up to two additional partons are used. 
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FIG. 17: The px distribution of the third hardest jet in inclusive Z/ 7 * production at the LHC as obtained by MCSNLD (red), 
PYTHIA (blue) and SHERPA (black). The dashed (dotted) curve corresponds to the prediction of SHERPA when matrix elements 
for up to two (three) additional partons are used. 


predicts a three jet rate of 9.6 pb, the samples with two (three) extra partons predict 16.3 (21.1) pb. However, this is 
not surprising keeping in mind that the LHC provides enough phase space to produce massive bosons in association 
with a multitude of high energetic jets that are best described by the corresponding matrix elements and that can 
not be appropriately described by parton shower emissions. 

To summarize: the predictions of MC0NLO and SHERPA agree fairly well for the shape of the boson transverse momentum 
and pseudo-rapidity distribution. Here, MCQNLO is of course superior in predicting the rate of inclusive Z and W 
production since it considers the corresponding production process at NLO in the coupling constant. This situation 
changes when studying the jets that potentially accompany the boson. As soon as more than one extra jet is considered 
SHERPA predicts significantly larger jet production rates and jet transverse momentum distributions that feature an 
enhanced population of the high-px region. Concerning PYTHIA it has to be stated that the shape of the boson 
transverse momentum and the boson pseudo-rapidity distribution differ significantly from the two other programs. 
This is directly related to the smaller amount of hard radiation produced by PYTHIA, clearly observed in the jet p± 
spectra. 
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V. SUMMARY 

In this publication, the previous validation of the merging procedure of matrix elements and the parton shower, as 
implemented in SHERPA, has been continued. Again, processes of the type pp ^ V + X, where V = lU^, have 
been chosen; this time, however, the analysis focused on the case of the CERN LHC rather than on the Fermilab 
Tevatron. Again, the merging procedure turned out to yield sufficiently stable results over a wide range of internal 
parameters, rendering it a predictive way of incorporating the full information available in tree-level matrix elements 
into multi-purpose event generators, as anticipated. In addition, when comparing the results obtained through SHERPA 
with those of a full next-to-leading order calculation, it again turned out that the results of SHERPA reproduce the 
essential features in the NLO shapes. However, it should be stressed that in SHERPA the total rates are still at leading 
order accuracy only. Nevertheless, the fact that SHERPA seems to reproduce the NLO shapes of the observables, the 
NLO rates can be recovered by simply multiplying with a constant AT-factor. When comparing the results of SHERPA 
with those of other event generators, some differences appear. Especially for observables sensitive to the correct 
treatment of multi-particle final states these differences have become significant, ranging up to orders of magnitude. 
In this study, SHERPA again proved its versatility in simulating high-multiplicity final states at collider experiments. 
Due to the merging procedure implemented in it, it provides a unique tool for the simulation of final states, where 
the proper treatment of the event topology is of great importance. 
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APPENDIX A: INPUT PARAMETERS AND PHASE-SPACE CUTS 

For all analyses, the PDF set cteqGl has been used, and as has been chosen according to the corresponding 
value of this PDF, namely as = 0.II8. For the running of the strong coupling constant, the corresponding two-loop 
equation has been employed. Jets or initial partons are restricted to the light flavor sector, namely g, u, d, s, and c. 
All flavors are taken to be massless. 


1. SM input parameters 


The SM parameters are given in the scheme: 

mw = 80.419 GeV, Tw = 2.06 GeV, 
mz = 91.188 GeV, Tz = 2.49 GeV, 

= 1.16639 X lO-*^ GeV“^ 

* 2/1 1 2/2 
sin = 1 - 

Os = 0.118. (AI) 


The electromagnetic coupling is derived from the Fermi constant according to 




sin 

TT 


(A2) 


The constant widths of the electroweak gauge bosons are introduced through the fixed-width scheme. CKM mixing 
of the quark generations is neglected. 


2. Cuts and jet criteria 

For all jet analysis the Run H k± clustering algorithm defined in has been used. The additional parameter of 
this jet algorithm is a pseudo-cone size /?, whose value has been chosen to R = 0.4. In addition jets have to fulfil the 
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following cuts on transverse momentum and pseudo-rapidity, 

pi^*>20GeV, < 4.5. 


(A3) 


For charged leptons the cuts applied are: 

> 15 GeV, < 2.4, mu > 15 GeV. (A4) 

No cut on missing transverse momentum has been applied. 

The final selection criteria correspond to the separation of the leptons amongst each other and with respect to the 
jets. 


ARij > 0.4 , ARu > 0.2. (A5) 

APPENDIX B: SETUPS FOR MC®NLO AND PYTHIA 

• The MC0NLO setup: The program version used is MCSNLD 2.31. The process number corresponding to pp —> 
Z/7* ^ e+e” -I- X production is IPR0C=-11351, for pp W~^ —> + X this equals IPR0C=-11461. In both 

cases consequently the underlying event has been switched off. The lepton pair in Zj^* production has been 
generated in a mass window of 


mz — SOr^ < TOee < TOz "P SOT^ (Bl) 

in the case of production the lepton-neutrino pair fulfils 

miY ~ SOTn/ < mei/^ ^ mw + SOTw (B2) 

The PDF set used is cteqOl. All other physics parameters that specify a run for MC0NLO have been left unchanged 
with respect to their default values. 

• The PYTHIA setup: The PYTHIA version used is 6.214. The process pp ^ Z/'y* +X is selected via the parameter 
MSUB(1)=1. The decay mode Zj^* e+e“ is picked by the settings MDME(182,1)=1 and MDME(170,1) = 1, all 
other decay channels have been disabled. The process pp -p A is turned on via MSUB(2) = 1. The decay 

mode IF+ —> e+t'e is chosen by MDME(206,1) = 1. It has proven to be convenient to increase the standard value of 
the shower starting scale in PYTHIA to = 14 TeV in order to produce a reasonable amount of high energetic 
QGD radiation. The parameter responsible for controlling the shower start scales is MSTP(68). For the above 
choice it has set to set to MSTP(68)=2. All other parameters specific for PYTHIA have been left unaltered, except 
that PYTHIA’s underlying event has been switched off. 
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